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ABSTRACT 

We present the results of a systematic numerical study of an alternative pro- 
genitor scenario to produce type la supernova explosions, which is not restricted 
to the ignition of a CO white dwarf near the Chandrasekhar mass. In this sce- 
nario, a shock-triggered thermonuclear explosion ensues from the collision of two 
white dwarfs. Consistent modeling of the gas dynamics together with nuclear 
reactions using both a smoothed particle and a grid-based hydrodynamics code 
are performed to study the viability of this alternative progenitor channel. We 
find that shock-triggered ignition and the synthesis of Ni are in fact a natural 
outcome for moderately massive white dwarf pairs colliding close to head-on. 
We use a multi-dimensional radiative transfer code to calculate the emergent 
broadband light curves and spectral time-series of these events. The synthetic 
spectra and lightcurves compare well with those of normal type la supernovae 
over a similar B-band decline rate and are broadly consistent with the Phillips 
relation, although a mild dependence on viewing angle is observed due to the 
asymmetry of the ejected debris. While event rates within galactic centers and 
globular clusters are found to be much too low to explain the bulk of the type 
la supernovae population, they may be frequent enough to make as much as 
a one percent contribution to the overall rate. Although these rate estimates 
are still subject to substantial uncertainties, they do suggest that dense stel- 
lar systems should provide upcoming supernova surveys with hundreds of such 
collision-induced thermonuclear explosions per year. 
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Introduction 



Type la supernovae (SNe la) are of major astrophysical relevance. They have acquired 
particular cosmological significance as a probe of the scale and geometry of the universe, 



providing the first evidence for its acceleration ( |Riess et al. 1998; Perlmutter et al. 1999 



Astier et al. 2006). These results depend crucially on the assumption that SNe la are 
standard candles. This assumption could be tested if the origins of SNe la are recognized. 
Knowledge of their nature is also of importance for understanding the metallicity evolution 
and star-formation history of galaxies. Yet, despite their relevance, no concensus on the 
nature of their progenitor systems has been reached. 

While there is broad agreement that the disintegration of a white dwarf (WD) in a 
thermonuclear explosion constitutes the supernova event itself, there are two main classes 
of competing models for the events that lead to the explosion. In the single-degenerate 



scenario, the exploding WD accretes from a non-degenerate stellar companion (Whelan & 



Iben 


1973; 


Nomoto 


1982) 



the double-degenerate scenario, the donor star is also a WD. The most commonly discussed 



progenitor system involves the coalescence of two CO WDs (Iben & Tut ukov 1984; Webbink 



1984), which after explosion should leave no remnant. There has been no conclusive proof 
to date that either scenario can lead to normal SNe la, nor has the evidence that the SN la 
rate is different for different stellar populations led to any firm conclusions. Therefore, any 
new observational or theoretical constraint on the progenitor systems is of great value. 

Here we present an alternative evolutionary scenario to produce a SNe la, which is 
not restricted to mass transfer in gravitationally bound double stellar systems. This new 
paradigm considers white dwarfs that reside within dense stellar systems where the stars are 
sufficiently close to each other to make collisions quite likely. The resultant shock compression 
could then lead to densities which exceed the threshold for pycnonuclear reactions so that 
thermonuclear runaway ensues. Understanding the feasibility of this channel for producing 
successful thermonuclear explosions as well as exploring the observational manifestations of 
such phenomena are the main purpose of this Letter. The layout is as follows. A concise 
summary of the numerical methods and the initial models is given in § 2. We describe 
the detailed hydrodynamic simulations in § 3, while the resulting broadband lightcurves 
and spectra together with a discussion of the relevance of this new progenitor channel to 
upcoming supernova surveys are presented in § 4. 
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Numerical Schemes and Initial Models 



As two stars approach each other to within a few stellar radii, their mutual gravitational 
interactions lead to the development of large-scale tidal distortions that substantially alter 
their global structures. If the trajectories of the stars bring them so close to each other that 
they experience a collision, the response of the stellar material to the impact is critical in 
understanding the future evolution of the system. As such, it is no longer appropriate to 
treat the stars as point masses, and a hydrodynamical description of the encounter becomes 
necessary. The outcome of a collision between two white dwarfs depends in an essential 
way on several factors: their masses and nuclear compositions, their relative speed, and the 
distance of closest approach. 

To study this problem, we use two complementary approaches: the Eulerian, adaptive- 



mesh hydrodynamics code FLASH (Fryxell et al. 2000), and a Lagrangian hydrodynamics 



code QRosswog et al. 2008; Rosswog et al. 2009) based on the Smoothed Particle Hydrody- 
namics (SPH) method ( Benz|[l990 ; Monaghan||2005 ; Rosswog||2009 ) . Both codes incorporate 
the Helmholtz equation of state ( Timmes fc Swesty||2000 ) and similar, small nuclear reaction 
networks tuned to correctly reproduce nuclear energy release ( Hix et aL]|1998 ; Timmes||l999 



Timmes et al.pOOO ), but they differ in their treatment of hydrodynamics and gravity. Using 
the same stellar models and impact conditions, this approach not only provides a classic 
code verification, but in particular allows to gain unique insight into the physics of white 
dwarf encounters. 

Some of the questions at the forefront of our attention are the effects of the initial 
nuclear composition and masses of the white dwarfs as well as the impact conditions. We 
have performed a large set of three-dimensional calculations. A detailed account of all models 
will be given elsewhere, here our focus is on central collisions. A summary of the performed 
calculations is given in Table [T] The stars with 0.4 M are instantiated as pure He, those 
with larger masses as a homogeneous mixture of 50% C and 50% O. All stars are initially 
cold (10 4 K), placed at a mutual distance of three times their combined unperturbed stellar 
radii and with the relative free-fall velocities of the corresponding point mass values. 



3. Shock- Triggered Thermonuclear Explosions from White Dwarf Impacts 

The relative velocity at contact is entirely dominated by the mutual gravitational at- 
traction, i.e. it is much larger than typical GC velocity dispersions, ctqc ~ 10 km/s, 
v rel = 4000 km/s (M t ot/1.2 M ) 1 / 2 (2 x 10 9 cm/(i?! + R 2 )) 1/2 > c s . The sound velocity c s in 
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the core of a 0.6 M WD is about 2600 km/s, and thus shocks are a natural result of a WD 
collision. Figure [T] shows the thermodynamic evolution of the most common combination of 
masses, 2 x 0.6 M . The snapshots of density and temperature in the upper two rows were 
obtained with SPH, the lower two are the result of FLASH (19-isotope network, minimum 
linear resolution 5 x 10 6 cm). The details of the collision differ slightly in the two simulation 
environments, which is evident in the shock geometries. However, the overall behavior is 
similar: shortly after contact, a discus-shaped, shock-heated region forms in which nuclear 
processing occurs. Since the ignition site is not coincident with the original central density 
peak of either WD, the shocks more quickly propagate through the shallow density gradient 
that is perpendicular to the direction of infall. As a result, the hot processed material first 
breaks out through a ring which lies in a plane that is parallel to the collisional plane. It is 
only when the rear of the stars have passed through the shock fronts that a significant overall 
expansion can set in. Apart from hydrodynamics and gravity, the codes also differ in the used 
reaction networks, the SPH code is coupled to a 7-species network ( Hix et aL]|1998 ) while the 



FLASH run uses a 19-isotope network |Timmes (1999). As a test of the energy generation 



accuracy we have used both networks along 1000 thermodynamic trajectories extracted from 
the SPH simulation. Maximum deviations in the resulting energy generation were 15%, while 
95% of the trajectories agreed to better than 5%. The mass fractions reported in this paper 
are all post-processed or direct results of the 19-isotope network. While the networks could 
be partly responsible for the shock structure differences the latter may also be influenced 
by the local resolution and details how burning directly in the shock is suppressed in both 
codes. A further difference between both runs is the spike-like feature at the rear side of 
the star in the FLASH simulation, this is an artifact of the rapid advection across the grid. 
Despite these differences our main conclusions are robust: the shocks trigger an explosion 
producing a substantial amount of 56 Ni (0.32 M for SPH, 0.16 M for FLASH). 

Mass-segregated environments may favour encounters of more massive WDs, though. 
We therefore show in Fig. [2] the outcome of a 2 x 0.9 M collision (run G; left: density and 
temperature, right: mass fractions). This event yields 0.66 M of 56 Ni, comparable to a 
typical SNIa. 



The topic of white dwarf collisions has been pioneered by |Benz et al. (1989). In their 



work they modeled the white dwarfs with an equation of state containing contributions 
from degenerate, relativistic electrons, from a Boltzmann gas of nuclei and from photons. 
Moreover, they coupled their smooth particle hydrodynamics code to a 14-isotope network. 
At that time their calculations were restricted to 5000 SPH particles resulting in a moderate 
numerical resolution. Although in some cases substantial nuclear burning took place, none of 
their models resulted in a complete disintegration of the white dwarf pairs. For comparison, 
we performed a test run of a headon collision between two 0.6 M WDs using 5000 SPH 
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particles in total, similar to their run 1. Consistent with their work, we find a surviving 
remnant and only 0.30 M of expelled material (0.09 M in their work). The degraded 
numerical resolution results in a one order of magnitude reduction of the released nuclear 
energy compared to our run D. The remaining differences between our test run and the 
results of Benz et al. are mainly due to the different equation of state, but they may to some 
extent also reflect the differences in the networks and the advances in the SPH method. Our 
overall results are similar to those of Raskin et al. (2009) which were submitted while our 
paper was under review. 



4. Discussion 

4.1. Lightcurves and Spectra 

Some white dwarf collisions should produce luminous light curves powered by the decay 
of radioactive 56 Ni synthesized in the explosion. To predict the observable emission, we post- 
processed select models using the radiative transfer code SEDONA ( Kasen et al.||2006 ). This 



code is 3-dimensional, time-dependent, multi-wavelength and includes a detailed treatment 
of the physics of 56 Ni decay and bound-bound line opacity. As initial conditions, we used 
the SPH results with post-processed abundances for the 2 x 0.6 M (run D) and 2 x 0.9 M 
(run G) collisions, at a simulation time late enough that the ejected material had reached 
the free, homologous phase of expansion. Given the axial symmetry of head-on encounters, 
the results where azimuthally averaged onto a cylindrical grid. 



In Figure [3j (right panel) we show synthetic spectra of the models, computed at the 
peak of the light curve (<~ 20 days after collision). The model spectra closely resemble 
those of normal Type la supernovae, with broad P-Cygni line features due to Si II, S II, 
and Ca II. This outcome is not surprising given that the ejecta compositional stratification 
(Figure [T]) is very similar to that of standard SNe la models, with an outer layer of interme- 
diate mass elements and an inner core of iron group elements. The asymmetry of the ejected 
debris introduces some variation of the spectrum with viewing angle, most prominently in 
the ultraviolet where the radiation transport is most sensitive to line blanketing opacity. 
The velocity of the supernova photosphere, as measured form the Doppler shift of the line 
absorption minima, is 13000 — 16000 km s _1 , similar to, though slightly higher than that 
observed in average SNe la. 

The light curves of the white dwarf collisions (Figure [3j center panel) also resemble 
those of Type la supernovae. As expected, the models which produced more 56 Ni have 
brighter light curves which decline more slowly. The calculated peak magnitudes and B- 
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band decline rates lie within the range of typical SNe la, and are consistent with the slope 
and normalization of the observed Phillips relation (Figure [3| left panel). This result is not 
totally unexpected, given that the principle parameter underlying the Phillips relation is the 



56 



Ni mass, which influences both the supernova luminosity and the ejecta opacity (Kasen 



& Woosley 2007, and references therein). On the other hand, the light curve width is also 



sensitive to the total ejected mass. In collision models, unlike most standard SN la scenarios, 
this value is not constrained to be the Chandrasekhar mass. Thus, although the particular 
models studied here roughly obey the Phillips relation, in detail white dwarf collisions could 
show small but systematic deviations. 



4.2. Diversity 

The explosion mechanism reported here is not tied to a particular mass scale and there- 
fore allows for considerably more diversity. As mentioned above, collisions between white 
dwarfs provide a pathway to ignite CO white dwarfs that completely disintegrate the WD 
pair. In contrast, for low-mass collisions (run A) or for collisions between a CO and a He 
WD a remnant remains, which, in the latter case, produces an identifiable outcome: a hot, 
high-speed (~ 1000 km/s) CO WD engulfed by a cloud of intermediate mass elements. 

The mechanism discussed here is found to work for the collision of two 0.4 M He WDs, 
but does not lead to an explosion in the case of 2 x 0.5 M . This is mainly the result of 
the different available fuel - helium burning releases more energy on the way to iron group 
elements (ene^Fe = 1-73 MeV/nucleon, eco^Fe = 0-98 MeV/nucleon). Due to the WD mass 
function, collisions between 0.6 M WDs are expected to dominate unless they occur in a 
strongly mass-segregated environments where more massive WDs would then be preferred. 
On average, however, this SN la channel will preferably ignite lighter WDs than standard 



channels (Hillebrandt & Niemeyer 2000; Podsiadlowski et al. 2008) and, as a result, the 



nucleosynthetic yields should be less neutron-rich due to the slower y e -evolution via electron 
captures. 



4.3. Detection Prospects 

In order to critically evaluate the outcome of these ideas, the frequency of such en- 
counters must be addressed. For a core-collapsed GC, the high densities in the core com- 
pletely dominate the collision rate. We assume the white dwarfs to be distributed homo- 
geneously within a spherical core of radius r c . We further assume that the total num- 
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ber density n w d and stellar velocity dispersion a c are constant within r c . Together with 
the dominating gravitational focusing, this means we can approximate the total collision 
rate as: is col = 20 Gyr" 1 /i/ 2 (n c /3 x 10 6 p C - 3 ) 2 (r c /0.1 pc) 3 (a c /10 kms" 1 )" 1 ([mi + 
777,2] /I M )(r co i/5 x 10 3 km), where fi < 1 is the fractional number of stars of type i 
within r c and we use the properties of the well-studied, proto-typical core- collapsed GC 



M15 (van den Bosch et al. 2006|. Fokker-Planck model fits to M15 predict the presence of 
a significant population of WDs with M > O.7M , so that fi > 0.5. By setting the distance 
of closest approach to the sum of the radii of the two WDs, r co \ = r\ + r 2 , an estimate of the 
collision rate i? co i « 10 2 / cc /i / 2 yr _1 Gpc ~ 3 , can be obtained by multiplying v co \ (per GC) 
with the average GC space density of n gc = 4.2 Mpc -3 (Brodie & Strader 2006). Here n. 



gc 



is derived by combining the number of GCs per galaxy with the galaxy luminosity density 
distribution and f cc is the fraction of core-collapsed clusters. This is most likely an under- 
estimation since, for example, the effect of binaries in GCs can increase v co \ by a moderate 
factor (~ 2, J. Fregeau, private communication). Moreover, if M 15 formed at higher initial 
concentration, it might have been in (or around) deepest core collapse for a longer time, 
significantly increasing the (average) v co \. Numerical experiments for the dominant 2 x 0.6 
M case indicate that about 20% of the collisions may result in explosions. 

Although these rates are still subject to significant uncertainties such as whether other 
GCs follow a similar core-collapse evolution, they indicate that white dwarf collisions in 
their dense cores are not unlikely and can contribute with a modest fraction to the SNe la 
population, whose event rates are estimated to be of order a few 10 4 yr _1 Gpc -3 (Cappellaro 



et al.|[l999 ). In addition, a number of collisions are also expected from ultra- compact dwarf 



galaxies (Hilker et al. 1999; Drinkwater et al. 2000, 2003), the hypercompact stellar systems 



that form when supermassive black holes are ejected from galactic centres by the action 



gravitational wave recoil (Merritt et al. 2009) and from more "typical" galactic centers. 



The transient sky at faint magnitudes is poorly known, but there are major efforts under 
way that would increase the discovery rate for type la supernovae from a few thousands to 
about hundreds of thousands per year. While the estimates given above are much too low to 
explain the bulk of the SNe la population, they may be frequent enough to provide upcoming 
supernova surveys with hundreds of collision-induced SNe la per year. 
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Table 1: Masses in M , initial separation a in Ri + R 2) densities and energies in cgs-units, 
"res." refers to the particle number for SPH calculations, and to maximum linear resolution 
in cm for FLASH. 



run 


masses 


a 


res. 


log(Pmax) 


2~max,9 


log(£ nuc ) 




remnant 












SPH 








A 


0.2, 0.2 


5 


2.0 x 10 5 


5.97 


2.5 


48.53 


0.044 


hot, He- WD 


B 


0.4, 0.4 


3 


2.0 x 10 6 


7.16 


4.0 


51.19 


0.80 


none 


C 


0.5, 0.5 


3 


1.0 x 10 6 


7.20 


4.8 


50.00 


0.21 


CO- WD in Ne-Mg-Si cloud 


D 


0.6, 0.6 


3 


2.0 x 10 6 


7.92 


8.9 


51.21 


1.20 


none 


E 


0.4, 0.9 


3 


2.5 x 10 6 


7.56 


3.4 


50.75 


0.40 


CO- WD in He-Si cloud 


F 


0.6, 0.9 


3 


2.0 x 10 6 


8.40 


7.9 


50.61 


0.30 


CO- WD in C-O-Si-Fe cloud 


G 


0.9, 0.9 


3 


1.0 x 10 6 


7.55 


6.3 


51.41 


1.80 


none 












FLASH 








H 


0.6, 0.6 


3 


4.9 x 10 6 


7.47 


5.5 


51.11 


1.20 


none 
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Fig. 1. — Comparison of density and temperature evolution of the central collision of two 0.6 
M CO white dwarfs. The upper two rows are the SPH-results, the lower two are produced 
by FLASH. The shown box length is 3 x 10 9 cm, limiting values colorbar (left to right): 
log(p) ma x = [7.12,7.05,6.91,6.84,6.70,6.64], T 9 , max = [1.49,4.79,3.99,3.65,3.39,3.24], 
log(p)min = 2 and T 9 ^ min = everywhere. 
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Fig. 2. — Collision between two 0.9 M white dwarfs: density and temperature (left), nuclear 
mass fractions (right). Note the different scales in both panels. 
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B-band decline rate AM 1 5 days since collision wavelength (Angstroms) 

Fig. 3. — Radiative transfer calculations of the light curves and spectra resulting from central 
collisions of 2 x 0.6 M (run D) and 2 x 0.9 M (run G) white dwarf pairs. The synthetic B- 
band light curves (central panel) closely resemble those of normal Type la supernovae. The 
peak brightness and decline rate of the light curves vary somewhat with the viewing angle 
(left panel), but are broadly consistent with the slope and spread of the observed Phillips 
relation (grey shaded band). The maximum light spectra (right panel) closely resemble that 
of the normal Type la supernova SN 1981B. 



